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Abstract 

Background: 5-Hydroxytryptamine receptor (5-HTR) and 5-hydroxytryptamine transporter (5-HTT) gene polymorphisms 
have been reported to be associated with susceptibility to obstructive sleep apnea syndrome (OSAS). The associations, 
derived from sporadic, inconsistent, small-sample-size studies, need to be evaluated further in a meta-analysis. 

Methods: Relevant studies were identified by searching PubMed, Embase, The Cochrane Library, China National Knowledge 
Infrastructure (CNKI), Wanfang, and Weipu. Eligible data were extracted from each included study. Odds ratios (ORs) were 
calculated using a fixed-effects or a random-effects model. The ORs and 95% confidence interval (CI) were used to assess the 
strength of the association between serotonergic gene polymorphisms and OSAS in the dominant and recessive models, as 
well as alleles. The Q statistic was used to evaluate homogeneity and Begg's test was used to assess publication bias. 

Results: Eight studies were finally included in the meta-analysis of the association between 5-HTR2A gene variants 
(including 1 02T/C and 1 438G/A), 5-HTT gene polymorphisms (including 5-HTT gene-linked promoter region (5-HTTLRP), and 
serotonin transporter intron 2 variable number tandem repeat (STin2VNTR) and OSAS risk. The G allele of 5-HTR2A 1438G/A, 
long 5-HTTLPR, and 10-tandem-repeats STin2VNTR were shown to increase OSAS susceptibility, with ORs of 2.33 (A vs. G, 
95% CI 1.48-3.66), 1.24 (L vs. S, 95% CI: 1.04-1.49), and 2.87 (10 vs. 12, 95% CI: 1.38-5.97), respectively. These significant 
differences were determined in both dominant and recessive models. Of the 5-HTR2A 1438G/A gene polymorphism, the AA 
genotype increased the OSAS risk, with an OR of 4.21 (95% CI: 2.83-6.25) in a recessive model in male OSAS patients, but no 
significant association was found in females. 

Conclusions: Our meta-analysis demonstrated that polymorphisms in the 5-HTR2A 1438G/A and 5-HTT genes contributed 
to susceptibility to OSAS. The A allele of the 1438G/A gene polymorphism is predominantly distributed in males and 
increased the OSAS risk significantly. 
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Introduction 

Because of its high morbidity and mortality, obstructive sleep 
apnea syndrome (OSAS) has recently become an important public 
health issue. It affects about 2% of females and 4% of males in 
middle-aged groups in Western counties [1]. Repeated partial or 
complete collapse of the pharynx during sleep, which results in 
apnea or hypopnea associated with oxygen desaturation and 
arousal from sleep, is the most characteristic feature of OSAS [2]. 
The pathophysiology of OSAS is associated with activation of a 
number of neural, humoral, thrombotic, metabolic, and inflam- 
matory disease mechanisms [3] . 

Serotonergic neurons innervating motoneurons increase their 
firing rates in response to respiratory challenges, and long-term 
facilitation of respiratory activity in response to hypoxia is 5- 



hydroxytryptamine (5-HT, serotonin)-dependent [4]. Thus, it is 
reasonable that the receptor and transporter of 5-HT are also 
involved in the pathophysiological process of OSAS. 

There is a large family of 5-HT receptors [5] . Among them, 5- 
HT receptor (5-HTR) 2A/2C has been most studied. The 5-HTR 
2A gene is located on chromosome 1 3ql4— q2 1 and the 5-HTR 2C 
gene is located in the q24 region of chromosome X [6-7]. Binding 
with the 2A/2C receptors, 5-HT plays a vital role in preventing 
glossocoma and sustaining patency of the upper airway [8]. A 
series of case control studies reported an association between 5- 
HTR 2A/2C gene polymorphisms and OSAS [9-14]. In 5-HTR 
2A, two polymorphisms have been reported: a 1438A/G 
polymorphism in the promoter region and an Mspl polymorphic 
site at position 102T/C [12]. A 796C/G polymorphism locus in 
the 5-HTR 2C gene has also been identified [13]. 
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The 5-HT transporter (5-HTT, serotonin transporter) gene is 
located in the 17ql 1 . 1 — 1 7ql 2 region of chromosome 17 [15]. Two 
gene polymorphisms, the 5-HTT gene-linked polymorphic region 
(5-HTTLPR) and the intron 2 variable numbers of tandem repeats 
(STin2VNTR), have been studied in connection with OSAS 
susceptibility because STin2VNTR can affect enhancer function 
and transcription of the 5-HTT gene and 5-HTTLPR can 
influence 5-HT reuptake [16-18]. 

Although the various studies have provided important seroto- 
nergic genetic information regarding susceptibility to OSAS [9- 
14,19,20], because of the limitations of small sample sizes in each 
study, low statistical power, and inconsistent results, they are 
incapable of demonstrating a reliable association. Also, the 
associations hinted at from sporadic cases and populations of 
differing ethnic origins need to be re-analyzed comprehensively. In 
this report, to further evaluate the associations between 5-HT 
receptor and transporter gene polymorphisms and the risk of 
OSAS, we performed a meta-analysis by combining the genotyp- 
ing data from eligible studies published to date. 

Materials and Methods 

This study was strictly followed the guidelines of the PRISMA 
(Preferred Reporting Items for Systematic Reviews and Meta- 
Analysis). 

Publication Search 

To identify all eligible studies that address the association 
between 5-HT receptor and transporter gene polymorphisms and 
OSAS, we conducted key word searches in the biomedical 
electronic databases, including PubMed, Embase, The Cochrane 
Library, China National Knowledge Infrastructure (CNKI), 
Wanfang, and Weipu: (obstructive sleep apnea hypopnea 
syndrome or OSAHS or obstructive sleep apnea syndrome or 
OSAS or obstructive sleep apnea or OSA or sleep apnea or apnea) 
and (5-hydroxytryptamine receptor or 5-HTR or serotonin 
receptor or 5-hydroxytryptamine transporter or 5-HTT or 
serotonin transporter) in combination with (gene or polymorphism 
or variants or alleles) were used as research terms. The final search 
was updated in June 2013. The search was limited to studies 
published in English and Chinese. 

Inclusion and Exclusion Criteria 

Case-control studies that met the following criteria were 
included: 1) published in peer-reviewed journals, 2) contained 
independent data, 3) had enough data to evaluate the relationship 
between gene polymorphisms and OSAS susceptibility, and 4) had 
enough data to determine the genotype distributions in both the 
case and control groups. The authors were contacted if there were 
queries regarding the studies. 

Studies were excluded for the following reasons: 1) lack of a 
suitable control group, 2) no usable data reported, 3) data were 
contradictory, 4) repeated data, and 5) abstracts or reviews. 

Data Extraction 

Two investigators (Drs. Xu and Guan) performed the data 
search, data screening, and data extraction independently. 
Disagreements regarding data screening and extraction were 
resolved by discussion. If there was still disagreement, a third 
reviewer (Prof. Yin) participated to resolve the issue. The following 
data were extracted from the original studies: first author, year of 
publication, country, sample size, age, gender, body mass index 
(BMI), apnea-hypopnea index (AHI), and genotype frequency. 



Statistical Methods 

The Review Manager software (ver. 5.0, Cochrane Collabora- 
tion, Oxford, UK), the SAS software (ver. 9.1.3, SAS Institute Inc., 
Gary, NC), and Stata (ver. 1 1 .0, Stata Corporation, College 
Station, TX) were used for the meta-analysis and statistical 
analyses. For genetic frequency distributions, Hardy-Weinberg 
equilibrium (HWE) was calculated with the SAS software to assess 
goodness of fit. Homogeneity was tested using Cochran chi- 
squared tests. The Mantel-Haenszel fixed-effects model or the 
random-effects model was used for the meta-analysis in the 
comparison of serotonergic gene polymorphisms [21—22]. Begg's 
test was used to evaluate publication bias. 

Results 

Search Results 

In total, 52 studies were identified by electronic and manual 
searching. At the stage of screening by the inclusion criteria, 38 
studies were excluded for the following reasons: animal studies, cell 
studies, reviews, or not correlated with OSAS. At the stage of 
screening by the exclusion criteria, after reading the full articles, 
six studies were excluded for the following reasons: two articles 
were repeated studies [23—24], one was an abstract [25], two 
articles lacked complete data [26-27], and one article had 
contradictory data [28]. Finally, eight studies met the exclusion 
criteria and were included in our systemic review and meta- 
analysis (Figure 1). Of the eight studies, six had data for 5-HTR2A 
102 T/C, five for 5-HTR2A 1438 G/A, three for 5-HTTLRP, 
and three for STin2VNTR (Figure 1, Table SI). 

Study Characteristics 

All eligible studies were published between 2005 and 2013. Five 
studies were from China [10,13,14,19,20], one from Turkey [9], 
one from Brazil [11], and one from Japan [12] (Table SI). The 
studied populations were mainly middle-aged and slightly over- 
weight in cases and controls, although one study lacked this 
information (Table SI). One study did not provide gender data 
[12]; the male proportion in the remaining populations ranged 
from 40.0% to 88.6% (Table SI). The BMI ranged from 
21.3±7.5 kg/m 2 to 29.3±3.7 kg/m 2 , although one study lacked 
these data [9]. The AHI of OSAS populations ranged from >5 to 
53. 9± 16.4 times per h; the AHI of control groups were all below 5 
(Table SI). All polymorphisms were identified by polymerase 
chain reaction (PGR) or PCR-restriction fragment length poly- 
morphism (PCR-RFLP) techniques. 

The allele frequencies of all OSAS and control samples were in 
HWE except one study (Table 1). The combined samples for 5- 
HT2A 102C/T [9-14], 5-HT2A 1438G/A [9-11,13,14], 5- 
HTTLPR [10,19,20], and STin2VNTR [10,19,20] polymor- 
phisms for the meta-analysis in case and control groups were 566, 
466, 593, and 593 and 728, 551, 479, and 479, respectively (Table 
SI). 

Meta-analysis of the Association between 5-HTR and 5- 
HTT Gene Polymorphisms and OSAS 

The association between gene polymorphisms and susceptibility 
to OSAS was evaluated in the total case and control populations 
combined from eligible studies of each gene locus. All genetic loci 
were tested in a dominant model, a recessive model, and a 
susceptibility allele versus insusceptibility allele model. 

Of the 5-HTR2A 102T/C polymorphisms, six studies contain- 
ing 566 case subjects and 728 controls were included in this meta- 
analysis [9-14]. The I 2 of the dominant models (TT+TC vs. CC), 
recessive model (TT vs. TC+CC) and T allele vs. C allele were 0 
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52 studies obtained from initial search 



38 studies were excluded by 
inclusion criteria 



14 studies were retained 



6 studies were excluded by 
exclusion criteria 



8 studies met all criteria 



6 studies on 
102T/D 9 - 14 ] 



5 studies on 

1438G/A[ 9 - 11 < 1 ^ 14 ] 



3 studies on 5- 
HTTLPR[ 1 °.19,20] 



3 studies on 5- 
HTTVNTRl 10 ' 19 ' 20 ! 



Figure 1. Flow chart of literature search and study selection. 6, 5, 3 and 3 studies were eligibly included in the meta-analysis of 5-HTR2A 102 
T/C, 5-HTR2A 1438 G/A, HTTLRP and HTTVNTR, respectively. 
doi:1 0.1 371 /journal.pone.0086460.g001 



Table 1. Genotypes and alleles distribution in OSAS and control groups. 







Genotypes 










Alleles 








HWE 


Study 


OSAS 






Control 






OSAS 




Control 




102C/T 


TT 


TC 


cc 


TT 


TC 


cc 


T 


C 


T 


c 




Bayazit Y, 2006 [9] 


17 


21 


17 


28 


54 


20 


55 


55 


110 


94 


Yes 


Chen H, 2013 [10] 


34 


56 


31 


29 


54 


22 


124 


118 


112 


98 


Yes 


De Carvalho T, 2013 [11] 


23 


66 


11 


20 


69 


11 


112 


88 


109 


91 


No 


Sakai K, 2005 [12] 


47 


90 


40 


25 


46 


29 


184 


170 


96 


104 


Yes 


Yin G, 2012 [13] 


52 


107 


51 


30 


51 


24 


NM 


NM 


NM 


NM 


Yes 


Zhu J, 2007 [14] 


18 


30 


17 


15 


28 


11 


66 


64 


58 


50 


Yes 


1438G/A 


AA 


AG 


GG 


AA 


AG 


GG 


A 


G 


A 


G 




Bayazit Y, 2006 [9] 


24 


19 


12 


25 


50 


27 


67 


43 


100 


104 


Yes 


Chen H, 2013 [10] 


74 


30 


17 


28 


28 


49 


178 


64 


84 


126 


Yes 


De Carvalho T, 2013 [11] 


35 


61 


4 


19 


71 


10 


131 


69 


109 


91 


No 


Yin G, 2012 [13] 


85 


82 


43 


17 


67 


21 


252 


168 


101 


109 


Yes 


Zhu J, 2007 [14] 


41 


15 


9 


15 


14 


25 


97 


33 


44 


64 


Yes 


5-HTTLPR 


SS 


SL 


LL 


SS 


SL 


LL 


S 


L 


S 


L 




Chen H, 201 3 [10] 


59 


35 


27 


55 


34 


16 


158 


84 


146 


64 


Yes 


Yue W, 2005 [20] 


51 


33 


20 


78 


54 


18 


135 


73 


210 


90 


Yes 


Yue W, 2008 [19] 


114 


106 


34 


173 


131 


34 


334 


174 


477 


199 


Yes 


STin2VNTR 


10/10 


10/12 


12/12 


10/10 


10/12 


12/12 


10 


12 


10 


12 




Chen H, 2013 [10] 


8 


24 


89 


1 


10 


94 


64 


178 


12 


198 


Yes 


Yue W, 2005 [10] 


3 


20 


81 


1 


14 


135 


26 


182 


16 


284 


Yes 


Yue W, 2008 [19] 


6 


46 


202 


3 


41 


294 


58 


450 


47 


629 


Yes 



Abbreviation: OSAS, obstructive sleep apnea syndrome; NM, not mentioned; S/L allele, short/long allele; 5-HTTLPR, 5-HTT gene-linked polymorphic region; STin2VNTR, 
serotonin transporter intron 2 variable number tandem repeat; HWE, Hardy-Weinberg equilibrium. 
doi:1 0.1 371 /joumal.pone.0086460.t001 
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and the p values of the Q_test were higher than 0.1, suggesting no 
significant heterogeneity among the studies, so we used a fixed- 
effects model to analyze the data. As shown in Table 2, meta- 
analysis of the dominant model (TT +TC vs. CC), recessive model 
(TT vs. TC +CC), and T allele vs. C allele revealed that no 
genotype or allele was associated with OSAS risk (Table 2, 
Figure 2A, B, C). 

For the 5-HTR2A 1438G/A polymorphism, five studies had 
466 case subjects and 551 controls [9-11,13,14]. The I 2 of the 
dominant models (AA+AG vs. GG) and the A allele vs. G allele 



were over 50%, and the p values of the Q test were 0.47 and 
<0.01, respectively, suggesting strong heterogeneity among the 
studies. So, these data were analyzed using a random-effects 
model. The meta-analysis revealed that the A allele was associated 
with OSAS risk with an OR of 2.44 (95% CI: 1.1 1-5.38) in the 
dominant model (AA+AG vs. GG). The A allele increased the 
OSAS risk with an OR of 2.33 (95% CI: 1.48-3.66) compared 
with the G allele (Table 2, Figure 2D, F). The I 2 of the recessive 
model (AA vs. AG+GG) was 0 and the p value of the Q test was 
less than 0.01, so we used a fixed-effects model to analyze the data. 



A. 1 02T/C, TT+TC vs. CC 



Bayazi V.20O6 


36 


58 


82 


102 


162% 


0 55 [0 26. I 16] 






Chan H.2013 


90 


121 


83 


105 


20 B'-, 


0.77 [0.41.1.431 






De Carvaino T.2013 




100 


89 


1 00 


9.0% 


1.00 [0.41, 2.43) 






Sakai K.2005 


137 


177 


71 


100 


18 8% 


1 40 [0.80. 2.44) 






YlnG.2012 


1» 


210 


81 


MB 


24.0% 


0 92 [0.53. 1.61) 






Zhu J.2007 


48 


66 


43 


54 


11 2% 


0 72 [0.30. 1.71) 






Total {95% CI) 




728 




566 


100.0% 


0.90 [0.69. 1.19) 


4 




Total events 
Hotorooonoitv Chi' = 


561 
4.67. Of = 


b(P = 


449 

.46): P ■ 


0% 











Test for overall effect: Z = 0 73 (P = 0.47) 



B. 102T/C, TT vs. TC+CC 



Odds R.il.o 



Favours OSAS Favours Control 



Otitis R.ilio 



Bayaz. Y.2006 


17 


BH 

60 


28 


102 


WfinM 
114% 


M-H Fi.r-d <tSKC.\ 

1.18 (0.58. 243) 


M-H Fli- 


rt 9S%CI 


Chen H.2013 


3d 


121 


29 


105 


19 1% 


1.02 [0.57. 1.84] 






De Carvalho r.2013 


23 


100 


20 


100 


nam 


1 19 [0 61. 2 35] 






Sakai K.2005 


47 


177 


25 


100 


20 1% 


1.08 [0.62. 1.90] 






Yin G.2012 


82 


210 


30 


105 


294% 


0.8210.49.1.391 






Zhu J.2007 


18 


65 


15 


54 


HON 


1.00 [0.44. 2.23] 






Total {95% CI) 




728 




S66 


100.0% 


1,02 (0.79. 1.32J 






Total events 

Molorooenciiv Ch ' 


191 
1.06. df- 


I* " 


147 

96V !' - 


0% 











Test tor overall effect: Z - 0.17 (P ■ 0.88) 



C. 102T/C, T allele vs. C allele 



Study or Sulwonn 

Bayazi Y.2006 


55 110 


110 


■■• 

20J 


■<■<■ 

18.4% 


M-H Fl*r,d 9S%CI 

045(0.54,1.38) 




1 W- CI 


Chen H.2013 


124 242 


112 


210 


24.9% 


0.92(0.64.1.331 






DeCarvalho T.2013 


112 200 


M 


2D0 


20.4% 


1 .06(0.72. 1.58] 






Yin G.2012 


184 354 


96 


2CO 


25.1% 


1.17(0.83. 1.66] 






Zhu J.2007 


66 130 


53 


106 


13 3v, 


0.89 (0.53. 148] 






Total (95% CI) 


1038 




922 


100.0% 


1.00 [0. S3. V20] 






Total events 

HBtsnvwiallu' f> = 


541 

1 74 rlf = & (P = 


485 

7B1 1' = 













Tesl for overall effect: Z = 0.03 {P = 0.98) 



D. 1438G/A, AA+AG vs. GG 



Sliitlv nrSnhnronn 

Bayazi Y.2006 


Fvpnta 

43 


95 


75 


102 


Wmohl 
20.4% 


M-H Random 95% CI 

1.29(0.59. 2.80] 


M-H. Rand 


nm 0S% CI 


Chen H.2013 


MM 


121 


56 


105 


21 6 : ; 


5.35 [2.82. 10.15) 






De Carvalho T.2013 


96 


100 


SO 


100 


16.1% 


2.67 (0.81.8.81) 






Yin G.2012 


:67 


210 


84 


105 


224% 


0.97 [0 54.1.74] 






Zhu J, 2007 


56 


65 


29 


54 


19.3% 


5.36(2.22, 12.98) 






Total (95% CI) 




851 




■J 66 


100.0% 


2.44 [1.11.5.38] 






Total events 


466 


i Ml I 


334 


p = n r 











Test for overall effect Z ■ 2.22 (P - 0 0: 



E. 1438G/A.AA vs. AG+GG 

OSAS Control Odds Ratio 

-Study flf.Sjjthjroup, faejittLJfltot-Etfoln,J^Uftflifllat BmfiajBtaa&CL. 



Odds Ratio 
M-H Flwrt 



Bayazi Y.20O6 


24 


55 


25 


102 


18.5% 


238 [1 19.4 79] 


Chen H.2013 


74 


121 


28 


105 


21.8% 


4.33 [2.46. 7.63] 


De Carvalho T.201 3 


35 


100 


19 


100 


23.1% 


2.30 [1.20.4.38] 


Yin G.2012 


85 


210 


17 


105 


25.3% 


3 52 [1 96.6.34) 


Zhu J, 2007 


41 


65 


15 


64 


11.3% 


4.44 [2 .04. 9.69] 


Total (95% CI) 




551 




■106 


100.0% 


3.31 [2.48. 4.40] 


Total events 


259 




104 








Heterogeneity: Chr* = 


3.53. df = 


4(P = 


)47); l» = 


OK 







Test for overall effect: Z ■ 8.19 (P <■ 0.00001) 



F. 1438G/A, A allele vs. G allele 



OSAS Control 
Study or Subgroup Events Total Evnnts Total Wnioht 



Odds Ratio 

M-H. Random. 85% CI 



Odds Ratio 
M-H. Random. 95% CI 



Bayazi Y.2006 
Chen H.2013 
DeCarvalho T.2013 
YlnG.2012 
Zhu J.2007 



97 130 



44 108 



19 3% 

20 6% 

20 5% 

21 6% 

ieo : : 



1.62(1.01.2.60] 
4.17(2.60. 6.21] 
1.59(1.06. 2.37) 
1.62(1.16. 2.26] 
4.28 (2.46. 7.42) 



Total (95% CI) 1102 
Total events 725 
Heterogeneity: Tau 1 = 0.22; Chi 1 = 22.96. df = 4 (P = 0.0001): 
Test for overall effect: Z - 3.68 (P ■ 0.0002) 



932 100.0% 2.33 [148. 3.68) 



G.LPRS/L, LL+SLvs. SS 



J. VNTR 10/12, 10/10+10/12 vs. 12/12 



Study or Subgroup 
Chen H.2013 


82 121 


50 


105 


■v. 

22.0% 


M-H Fur-d 95% CI 

1.16 [0.89.145] 


M-H Fh» 


|d 95% CI 


Yue W.2005 


53 104 


72 


150 


244% 


1 13(0 68. 1 86] 






Yue W.2008 


140 254 


165 


336 


53.6% 


1.29 [0.93. 1.78) 




m 


Total (95% CI) 


479 




593 


100.0% 


1.22 [0.96,1.55] 




♦ 


Total events 
Hotcroaene.1v. Chi' = 


255 

0.24.df>2(P> 


287 
.881: 1* « 


0'-, 











Sturiv nr Ri.hnmnn 


Fvpnts Tnt.il 






Weiohl 


M-H Fi>pd qS% CI 


M-H Fix 


rl «IS%CI 


Chen H.2013 


32 121 


11 


105 


18.0% 


3.07 [1,46. 6.46J 






Yue W.2005 


23 104 


15 


150 


19.8% 


2.56 [1.26. 5.18J 






Yue W.20O8 


52 254 


44 


338 


62.2% 


1.72 [1.11. 2.67] 






Total (95% Ci) 


479 




593 


100.0% 


2.13(1.33,2.96) 




♦ 


Total events 

Hiilfiri-viorviplu Ch, 1 - 


107 

7 in rtl ■ ■> tO m 


70 


Ret 











Test for overall effect Z = 1.60 (P a 0.11) 



Favours OSAS Favours control 



H. LPRS/L, LLvs.SL+SS 



Slufjy ei Sufaaiaua Events Total Events Total Weight M-H, Fiatd. 3Ki CI 



Odds Ratio 

v-i i y ;'„■ 



Chen H.2013 
Yue W.2005 
Yue W.2008 



34 254 



105 26 4% 1.60 [0 81. 3.16] 

150 23.6% 1.7510 97.349] 

338 50.1% 1.38 [0.83, 2.291 

593 100.0% 1.52 [1.07. 2.16] 



Total (95% CI) 479 
Total events 81 68 

Heterogeneity: Ch? = 0.31. df = 2 (P = 0,86): P = 0% 
Test for overall effect Z = 2.36 (P = 0.02) 

I. LPRS/L, S allele vs. L allele 

OSAS Control Odds Ratio 



Test for overall effect: Z - 4.47 (P < 0.00001 > 

K. VNTR 10/12, 10/10 vs. 10/12+12/12 

OSAS Control Odds Ratio 

Study or Subgroup Events Total Events Total Weight M-H. Fixed. 95% CI 



Odds Ratio 

M-H.Fhted.95%Ct 



Chen H.2013 
Yue W.2005 
Yue W.2008 



2V. 
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Test for overall effect: Z = 2.72 (P = 0.007) 
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L. VNTR 10/12, 10 allele vs. 12 allele 
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210 
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73 


208 


90 


300 
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Yue W.2008 
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1.25 [0.98. 1.60) 
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Test for overall effect Z = 2.83 (P = 0.005) 



Figure 2. Forest plots for the associations between 5-HTR2A 102 T/C, 5-HTR2A 1438 G/A, HTTLRP and HTTVNTR gene 
polymorphisms and OSAS risk. The squares correspond to the study specific OR and 95% CI. The area of the squares reflects the weight. The 
diamond represents the summary OR and 95% CI. 
doi:1 0.1 371 /journal. pone.0086460.g002 
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Table 2. Meta-analysis of the association between gene polymorphisms and OSAS risk. 







Genotypes 


No. 


Test of association 




Test of heterogeneity 




Polymorphism 






OR 


95% CI 


P value 


Model 


Q 




I 2 


102T/C 


TT+TC vs.CC 


6 


0.90 


0.69-1.19 


0.47 


F 


4.67 


0.46 


0 




TT vs. TC+CC 


6 


1.02 


0.79-1.32 


0.86 


F 


1.06 


0.96 


0 




T allele vs. C allele 


5 


1.00 


0.83-1.20 


0.98 


F 


1.74 


0.78 


0 


1438G/A 


AA+AG vs. GG 


5 


2.44 


1 .1 1 -5.38 


0.03 


R 


20.76 


0.47 


81 




AA vs. AG+GG 


5 


3.31 


2.48-4.40 


<0.01 


F 


3.53 


<0.01 


0 




A allele vs. G allele 


5 


2.33 


1 .48-3.66 


<0.01 


R 


22.96 


<0.01 


83 


5-HTTLPR LIS 


LL+SL vs. SS 


3 


1.22 


0.96-1.55 


0.11 


F 


0.24 


0.86 


0 




LL vs. SL+SS 


3 


1.52 


1.07-2.16 


0.02 


F 


0.31 


0.88 


0 




L allele vs. S allele 


3 


1.24 


1.04-1.49 


0.02 


F 


0.02 


0.99 


0 


STin2VNTR 10/12 


10/10+10/12 vs.12/12 


3 


2.13 


1.53-2.96 


<0.01 


F 


2.1 


0.35 


5 




10/10 vs.10/12+12/12 


3 


4.10 


1.48-11.35 


<0.01 


F 


0.65 


0.72 


0 




10 allele vs. 12 allele 


3 


2.87 


1.38-5.97 


<0.01 


R 


10.18 


<0.01 


80 



Abbreviation: L or S, long or short 5-HTT gene-linked polymorphic region;1 0 or 1 2, serotonin transporter intron 2 variable number tandem repeat; F, fixed-effects model; 
R, random-effects model; OR, odds ratios; 95%CI, 95% confidence interval. 
doi:1 0.1 371 /journal.pone.0086460.t002 



As shown in Table 2 and Figure 2E, the A allele was associated 
with OSAS risk with an OR of 3.31 (95% CI: 2.48-4.40) in the 
recessive model (AA vs. AG+GG). 

Of the long or short 5-HTTLPR, three studies were pooled for 
totals of 593 and 479 subjects in the OSAS and control groups 
[10,19,20], respectively. The I 2 of the dominant model (LL+SL vs. 
SS), recessive model (LL vs. SL+SS), and L allele vs. S allele were 0 
and p values of the Q_test were all more than 0.1. Meta-analysis 
with a fixed-effects model showed that the LL genotype increased 
the OSAS risk at an OR of 1.52 (95% CI: 1.07-2.16) when 
compared with a population with LS or SS genotypes. A single 
allele L increased the OSAS risk at an OR of 1.24 (95% CI: 1.04- 
1.49) compared with the allele S (Table 2, Figure 2G and I). The 
dominant model (LL+SL vs. SS) showed no significant association 
with OSAS risk (Table 2, Figure 2H). 

For the STin2VNTR gene polymorphism, three studies were 
pooled with 593 and 479 subjects in the OSAS and control groups, 
respectively [10,19,20]. The I 2 of the dominant model (10/10+10/ 
12 vs. 12/12) and recessive model (10/ 10 vs. 10/12+12/12), were 
5 and 0, and the p values of the 0_ test were 0.35 and 0.72, 
respectively. Meta-analysis using a fixed-effects model showed that 
10 tandem repeats increased the OSAS risk with an OR of 2.13 
(95% CI: 1.53-2.96) and 4.10 (95% CI: 1.48-1 1.35) in a dominant 
model (10/10+10/12 vs. 12/12) and recessive model (10/10 vs. 
10/12+12/12), respectively (Table 2, Figure 2J, K). The 10 
tandem repeats variant increased the OSAS risk at an OR of 2.87 
(95% CI: 1.38-5.97) compared with 12 tandem repeats (Table 2, 
Figure 2L). 

Genotype Distribution of 5-HTR 2A and 5-HTT Gene 
Polymorphisms Sorted by BMI, AHI and Gender 

To further evaluate the association between 5-HTR 2A and 5- 
HTT gene polymorphisms and OSAS susceptibility, we sorted the 
available data according to BMI, AHI, and gender. The genotype 
distributions appeared equal within groups categorized by BMI 
and AHI indexes, consistent with previous reports (data not 
shown). For gender, two [9,1 1], four [9-11,14], one [19], and one 
[19] studies addressed the 5-HTR2A 102C/T, 5-HTR2A 1438G/ 
A, 5-HTTLPR L/S, and STin2VNTR 10/12 polymorphisms. 



The genotype distribution of 5-HTR2A 102C/T, 5-HTTLPR L/ 
S, and STin2VNTR 10/12 polymorphisms showed no significant 
difference between males and females (Table 3), while the AA 
genotype of 5-HTR 2A 1438G/A in males was significantly more 
common than in females. The AA genotype increased the OSAS 
risk with an OR of 4.21 (95% CI: 2.83-6.25) in a recessive model 
in male OSAS patients; however no significant association was 
found in females (Figure 3, Table 3). 

Publication Bias 

Begg's test was used to evaluate publication bias in the 
literature. The shapes of the funnel plots revealed no obvious 
asymmetry (not shown). We then used Begg's test to assess 
statistical evidence of funnel plot symmetry, but no significant 
publication bias was indicated (Table 4). 

Discussion 

Although the genetic basis of OSAS is poorly understood, 
studies to date have focused on relationships between gene 
polymorphisms and susceptibility to OSAS [9-14,19,20,29]. 
Among genetic variants, polymorphisms in 5-HTR and 5-HTT 
have been investigated most frequently. Thus, the associations 
suggested from sporadic, inconsistent, small-sample-size studies 
needed to be further evaluated in a systematic review and meta- 
analysis. In this report, the combined number of samples for 5- 
HT2A 102C/T, 5-HT2A 1438G/A, 5HTTLPR, and 
STin2VNTR polymorphisms in cases and controls were 566, 
466, 593, and 593 and 728, 551, 479, and 479, respectively. Our 
results revealed that the A allele of 5-HT2A 1438G/A, the long 5- 
HTTLPR, and 1 0 tandem repeats of STin2 VNTR increased the 
OSAS risk and that 5-HT2A 102T/C was not associated with 
OSAS. All risk alleles displayed increased susceptibility to OSAS 
in both dominant and recessive models except the long 5- 
HTTLPR, which displayed no significant association with OSAS 
risk in a dominant model (LL+SL vs. SS); this may be explained by 
an insufficient sample size or unknown factors. Of the 5-HTR2A 
1438G/A polymorphism, the AA genotype increased the OSAS 
risk with an OR of 4.21 (95% CI: 2.83-6.25) in a recessive model 
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A. 1438G/A,AAvs.AG+GG(male) 



OSAS 



Control 



Odds Ratio 



Study or Subaroup 


Events 


Total 


Events 


Total 


Weiqht 


M-H. Fixed. 95% CI 


Bayazit Y.2006 


19 


33 


14 


55 


18.3% 


3.97 [1.59, 9.96] 


Chen H.2013 


65 


99 


25 


89 


37.1% 


4.89 [2.63, 9.11] 


De Carvalho T.2013 


33 


73 


9 


40 


26.1% 


2.84 [1.19,6.81] 


Zhu J.2007 


36 


54 


12 


42 


18.5% 


5.00 [2.08, 12.01] 


Total (95% CI) 




259 




226 


100.0% 


4.21 [2.83, 6.25] 


Total events 


153 




60 








Heterogeneity: Chi 2 = 


1.17, df = 


3(P = 


D.76); I 2 = 


0% 







Odds Ratio 
M-H. Fixed. 95% CI 



Test for overall effect: Z = 7.1 1 (P < 0.00001 ) 



0.01 0.1 1 10 100 

Favours OSAS Favours control 



B. 1438G/A.AA vs. AG+GG (female) 



OSAS 



Control 



Odds Ratio 



Study or Subgroup Events Total Events Total Weight M-H. Fixed. 95% CI 



Odds Ratio 
M-H. Fixed. 95% CI 



Bayazit Y.2006 


5 


22 


11 


47 


35.0% 


0.96 [0.29, 3.21] 


Chen H.2013 


9 


22 


4 


16 


17.7% 


2.08 [0.50, 8.55] 


De Carvalho T.2013 


2 


27 


10 


60 


37.2% 


0.40 [0.08, 1.97] 


Zhu J.2007 


5 


11 


3 


12 


10.1% 


2.50 [0.43, 14.61] 


Total (95% CI) 




82 




135 


100.0% 


1.11 [0.56, 2.19] 


Total events 


21 




28 








Heterogeneity: Chi 2 = 


3.20, df = 3 


(P = 0.36); I 2 = 


6% 






Test for overall effect 


Z = 0.29 (P 


= 0.77) 











0.01 0.1 1 10 100 

Favours OSAS Favours control 



Figure 3. Forest plots for the association between 5-HTR2A 1438 G/A and OSAS risk in male and female groups. The squares 
correspond to the study specific OR and 95% CI. The area of the squares reflects the weight. The diamond represents the summary OR and 95% CI. 
doi:1 0.1 371 /journal.pone.0086460.g003 



in male OSAS patients; however no significant association was 
found in the female population. Other gene polymorphisms, such 
as 5-HTR2A 1421C/T, 1678A/C, 24563A/C, 61472A/C, and 
5-HTR2C 796C/G, failed in the meta-analysis due to the 
limitations of eligible studies or insufficient data. We attempted 
to categorize the data according to BMI, AHI, and gender to 
further evaluate the association between gene polymorphisms and 



OSAS. However, this was not possible due to the limited quantity 
of data available. 

Although race is a major factor in the distribution of gene 
polymorphisms, and our study data were from Japan [12], China 
[10,13,14,19,20], Turkey [9], and Brazil [11], we found no 
distribution differences in the gene polymorphisms evaluated. We 
could not conduct a haplotype analysis in this study because we 



Table 3. Genotype distribution of 5-HTR 2A and 5-HTT gene polymorphisms in OSAS patients subgrouped by gender. 





Study 


Male 




Female 




Chi-square with 
Yates' correction 


Two-tailed p value 


102C/T 


CC 


TT+TC 


CC 


TT+TC 






Bayazit Y, 2006 [9] 


8 


25 


9 


13 






De Carvalho T, 2013 [11] 


9 


64 


2 


25 






Pooled 


17 (16%) 


89 (84%) 


11 (22.4%) 


38 (77.6%) 


0.548 


0.459 


1438G/A 


AA 


GG+GA 


AA 


GG+GA 






Bayazit Y, 2006 [9] 


19 


14 


5 


17 






Chen H.2013 [10] 


65 


34 


9 


13 






De Carvalho T, 2013 [11] 


33 


40 


2 


25 






Zhu J.2007 [14] 


36 


18 


5 


6 






Pooled 


153 (59.1%) 


106 (40.9%) 


21 (25.6%) 


61 (74.4%) 


25.586 


<0.0001 


5-HTTLPR LIS 


LL 


SS+SL 


LL 


SS+SL 






Yue W, 2008 [19] 


30 (13.6%) 


190 (86.4%) 


4 (11.8%) 


30 (88.2%) 


0.001 


0.978 


STin2VNTR 10/12 


10/10 


12/12+10/12 


10/10 


12/12+10/12 






Yue W, 2008 [19] 


4 (1.8%) 


216 (98.2%) 


2 (5.9%) 


32 (94.1%) 


0.715 


0.398 



Abbreviation: 5-HTTLPR US, L or S, long or short 5-HTT gene-linked polymorphic region; STin2VNTR 10/12, 10 or 12, serotonin transporter intron 2 variable number 
tandem repeat. The genotype distribution was expressed as frequency or frequency with percentage. 
doi:1 0.1 371 /joumal.pone.0086460.t003 
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Table 4. Tests for publication bias in overall population. 



Polymorphism 


Comparision 


Bgger's test (p 
value) 


102T/C 


TT+TC vs.CC 


0.71 




TT vs. TC+CC 


0.71 




T allele vs. C allele 


0.73 


1438G/A 


AA+AG vs. GG 


0.81 




AA vs. AG+GG 


0.81 




A allele vs. G allele 


0.81 


5-HTTLPR LIS 


LL+SL vs. SS 


1.00 




LL vs. SL+SS 


1.00 




L allele vs. S allele 


1.00 


STin2VNTR 10/12 


10/10+10/12 vs.12/12 


1.00 




10/10 vs.10/12+12/12 


1.00 




10/10 allele vs. 12/12 allele 


1.00 



Abbreviation: 5-HTTLPR US, long or short 5-HTT gene-linked polymorphic 
region; STin2VNTR 10/12, 10 or 12 tandem repeats of serotonin transporter 
intron 2 variable number tandem repeat. 
doi:1 0.1 371 /journal.pone.0086460.t004 



could not collect original haplotype data. For the 5-HTR2A 
1438G/A and STin2VNTR meta-analyses, heterogeneity existed 
among the studies included, which, together with the analysis of 
the genotype distribution of 5-HTR 2A and 5-HTT gene 
polymorphisms according to gender, suggests that the heteroge- 
neity may have been caused by gender imbalances in the different 
studies. 

During the preparation of this manuscript, a systematic review 
and meta-analysis of the association between 5-HT2A receptor 
polymorphisms and risk of OSAS was published [30]. The 
differences between this report and our study are as follows: our 
study included not only 5-HTR gene polymorphisms but also 
5HTTLPR and STin2VNTR gene polymorphisms, and our study 
first demonstrated that gender was an important cofactor for 
1438G/A allele distribution; the AA genotype increased the OSAS 
risk with an OR of 4.21 (95% CI: 2.83-6.25) in a recessive model 
in male OSAS patients, but no significant association was found in 
females. 
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